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CHEMICAL  INTERACTIONS  AND  LIGHT  EMISSIONS 
FROM  VENTED  SPECIES 


I.  INTRODUCTION 

Airglov  emissions  may  be  induced  by  the  interaction  of  atmospheric 
constituents  with  vapors  released  from  space  vehicles.  Such  emissions  can 
be  produced  by  a  number  of  chemical  mechanisms.  Previous  studies  have 
considered  neutral-neutral  energetic  reactions  as  the  causative  mechanism 
for  vibrational  or  electronic  excitation  of  molecular  species  [Fraser  et 
al.,  1988].  ‘'This  study  considers  neutral-plasma  reactions  as  veil  as 
neutral-neutral  reactions  as  the  sources  for  excited  species.  To  provide 
an  example  of  the  type  of  interactions  which  can  occur,  we  model  the 
chemical  reactions  and  subsequent  emissions  resulting  from  water  (H20)  and 
methyl  alcohol  (CH^OH)  outgassed  from  a  spacecraft  in  the  upper  atmosphere 
above  300  km  altitude.  The  plasma  interactions  involve  0+  and  H+  thermal 
ions  and  low  energy  (<0.2  eV)  electrons  of  ionospheric  origin.  Line  and 
continuum  emissions  extending  from  the  infrared  to  the  extreme  ultraviolet 
are  of  interest.  The  examples  in  this  study  are  restricted  to  releases  of 
a  few  kilograms.  The  techniques  described  here  can  be  extended  to 
determine  the  emissions  caused  by  the  release  of  materials  besides  1^0  and 

ch3oh. 


Manuscript  approved  August  7,  1989. 


II.  CHEMICAL  REACTIONS 

A  released  species  initially  reacts  with  0+  or  H+  ions  to  produce  a 
posi tive-ion-intermediary  (PII)  [Bernhardt,  1987].  The  PII's  are  destroyed 
by  a  dissociative  recombination  reaction  which  may  yield  excited  species 
which  subsequently  decay  by  the  emission  of  a  photon.  The  complete  process 
is  summarized  by  the  reactions: 

XY  +  Z+  -»  XZ+  +  Y  rate:  ~  10~9  cm3/s  (1) 

XZ+  +  e"  -»  X*  +  Z  rate:  k2  ~  10~7  cm3/s  (2) 

X*  -»  X  +  h\>  rate:  A  ~  10-3  to  10^  s~*  (3) 

where  XY  is  a  generic  symbol  for  the  released  chemical,  Z+  represents  0+  or 
H+,  XZ+  is  the  positive  ion  intermediary  and  X*  is  the  excited  species 

which  radiates  a  photon  hv. 

A.  Ion-Molecule  Reactions 

The  primary  ion-molecule  reactions  for  the  ionospheric  releases  are 


0+  +  H20  -»  H20+(v<6) 

+  0  +  1.02 

eV 

^1A  "  x 

(A) 

H+  +  H20  -»  H20+(v<6) 

+  H  +  1.00 

eV 

k1B  =  (2.4  x  10~9? ) 

(5) 

0+  +  CH30H  -*  CH30+  +  0H(v)  +  ?  eV  (70%) 

-*  CH30H+  +  0*  +  2.78  eV  (25%)  kJC  =  1.9  x  10"9  (6) 

->  CH20+  +  H20  +  6.90  eV  (  5%) 

H+  +  CH30H  ->  CH30+  +  H2 

-»  CH4+  +  0H(v<2)  +  1.5  eV  kJD  =  (1.9  x  10-9?)  (7) 

-»  CH30H+  +  H  +  2.76  eV 

3 

The  rate  constants  (in  cm  /s)  are  from  laboratory  studies  at  temperatures 
of  300  K  [Albritton,  1878;  Ikezoe  et  al.,  1987].  The  (0+,  H-O)  reaction 
rate  (k^)  is  nearly  constant  for  energies  up  to  5  eV  [Murad,  private 
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communication,  1988]. 

Secondary  ion-molecule  reactions  occur  between  the  reaction  products 
and  the  injected  or  background  species. 


H20+ 

+  H20  ->  H30+  +  0H(v<2)  +  1.17  eV 

klE  = 

1.67  x 

10"9 

(8) 

h20+ 

+  °2  ->  02+  +  H20  +  0.54  eV 

klF=  4 

o  in-10 

.3  x  10 

(9) 

0+  + 

OH  ■*  0H+(v<l)  +  0  +  0.45  eV 

klG  ~ 

1.6  x 

10~9 

(10) 

0+  + 

H2  -♦  0H+(v<l)  +  H  +  0.37  eV 

klH  = 

2.0  x 

10-9 

(11) 

0+  + 

CH20  -»  CH0+  +  0H(v<9)  +  4.49  eV  (40%) 
-*  CH20+  +  0*  +  2.74  eV  (60%) 

kll  ‘ 

3.5  x 

10"9 

(12) 

0+  + 

CH4  •*  CH4+  +  0  +  1.01  eV  (89%) 

-♦  CH3+  +  0H(v<9)  +  3.78  eV  (11%) 

klJ  = 

1.0  x 

10-9 

(13) 

0+  + 

CH3  -»  CH3+  +  0*  +  3.77  eV 

klK  * 

1.0  x 

10-9 

(14) 

o 

+ 

+ 

CH2  -»  CH2+  +  0*  +  3.22  eV 

klL  “ 

1.0  x 

10-9 

(15) 

0+  + 

CH  -♦  CH+  +  0*  +  2.49  eV 

klM  = 

1.0  x 

10'9 

(16) 

h30+ 

+  ch3oh  -»  ch5o+  +  h2o 

klN  = 

2.2  x 

10~9 

(17) 

ch3o 

f  +  h2o  -»  h3o+  +  CH20(v) 

k10  = 

2.3  x 

10~9 

(18) 

ch2o 

f  +  H20  H30+  +  CHO  +  0.8  eV 

klP  = 

2.6  x 

10"9 

(19) 

<V 

+  H20  -►  H30+  +  CH3  +  1.87  eV 

klQ  = 

2.6  x 

10-9 

(20) 

<V 

+  CH30H  ■+  CH40+  +  CH4  +  1.8  eV  (60%) 
CH50+  +  CH3  (40%) 

klR  = 

3.0  x 

10'9 

(21) 
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CH30H+  4  CH30H  4  CH50+  4  CH30 

kls  =  2.5  x  10  9 

(22) 

CH50+  4  CH30H  4  H+-CH3OCH3  4  H20 

k1T  .  1.0  X  10-10 

(23) 

3 

where  all  rate  constants  are  in  cm  /s.  For  all  ion-molecule  reactions 

which  may  yield  excited  oxygen  (0*),  it  is  assumed  that  102  is  in  the 

1  3 

metastable  0(  D)  state  and  that  902  is  in  the  ground  0(  P)  state.  The 

value  for  k^  in  (10)  is  an  estimate  based  on  similar  reactions. 

If  known,  the  excess  energy  of  each  reaction  is  indicated.  If  this 

energy  is  much  less  than  the  injection  kinetic  energy  relative  to  0+  or  H+, 

then  the  reaction  rates  may  be  significantly  different  than  the  thermal 

rates  shown.  For  injection  velocities  less  than  5  km/s,  the  kinetic  energy 

of  0+  is  less  than  2.08  eV  and  of  H+  is  less  than  0.13  eV.  All  of  the 

chemical  rate  constants  k1A  through  k^g  are  valid  for  these  energies. 

The  excess  energy  also  indicates  which  species  may  be  electronically 

excited.  The  energy  of  reaction  plus  the  orbital  kinetic  energy  may 

combine  to  leave  the  reaction  products  in  excited  states  not  attainable 

from  low  energy  releases. 

B.  Dissociative  Recombination  Reactions 

The  amount  of  energy  released  during  dissociative  recombination  is 
given  by  Eq  =  I(XZ)  -  D(XZ)  where  I  and  D  are  the  respective  ionization 
potential  and  dissociation  energy  of  molecule  XZ  in  reaction  (2).  The  ion- 
electron  reactions  for  the  modified  plasma  are 

H20+  4  e"  4  OH(X2n.,v<9)  4  H  +  7.45  eV  (552) 

4  OH(A2Z+,v<9)  4  H  4  3.40  eV 

->  0(X3P)  4  H2  4  7.64  eV 

4  0(A1D)  4  H2  4  5.67  eV  (102)  kJA  =  6.5  x  10~7  (24) 

4  0(B1S)  4  H2  4  3.00  eV 
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(35%) 


-*•  0(X3P)  +  2H  +  3.00  eV 
-+  0(A1D)  +  2H  +  1.03  eV 

H30+  +  e“  -»  H20(X1A1  f v)  +  H  +  6.29  eV  (14%) 

0H(X2nifv<9)  +  H2  +  5.63  eV  (74%)  k2B  =  6.3  x  10~7  (25) 

0B(A2E+,v<5)  +  H2  +  1.58  eV 

->  0B(X2n.,v<2)  +  2B  +  1.09  eV  (12%) 

02+  +  e"  ->  0*  +  0*  +  6.90  eV  k2c  =  1.9  x  10“7  (26) 

0H+  +  e~  ->  0*  +  H  +  8.74  eV  k2D  -  7.5  x  10~8  (27) 

CBO+  +  e'  ■»  CO*  +  H  +  7.72eV  k2£  =  1.1  x  10~7  (28) 

CB20+  +  e-  ->  CBO  +  B  +  7.11  eV  k2F  "  7  (29) 

CB30+  +  e"  •*  CB20(v)  +  B  k2G  “  7  (30) 


CH3OB+  +  e~  ■»  CB20(v)  +  H2  +  9.92  eV 

4  CH20(v)  +  H  +  H  +  5.4  eV  k2B  =  7  (31) 

-»  CH2  +  B20*  +  6.83  eV 


CH50+  +  e_  4  CB4  +  OB 

k2l 

X 

00 

00 

II 

h** 

O 

I 

(32) 

C2B?0+  +  e~  4  C2B60  +  B 

k2J 

_  ? 

(33) 

CH . +  +  e"  4  CH3  +  B  +8.16  eV  (50%) 

4  CB2  +  B2  +  7.90  eV  (50%) 

u: 

CM 

.X 

II 

00 

o 

X 

10"7 

(34) 

CB3+  +  e'  4  CB2  +  H  +  4.99  eV  (50%) 
-»  CB  +  fl2  +  5.08  eV  (50%) 


k2L  =  8.0  x  10'7  (35) 


CB2+  +  e-  4  CB  + 

H  +  6.03  eV 

k2M  =  6.0  x 

10  7 

(36) 

CH+  +  e-  4  c  +  H 

+  7.63  eV 

k2N  =  x 

J— » 

o 

1 

(37) 

5 


3  1/2 

where  the  rate  constants  (cm  /s)  all  scale  as  (300/T^)  ,  Tg  is  the 

electron  temperature  and  the  unknown  rate  constants  have  values  around  6  x 
7  3 

10  cm  /s  [Mitchell  and  McGowan,  1983].  The  percentages  for  the  products 

for  dissociative  recombination  of  H20+  and  H30+  are  taken  from  Rowe  et  al. 

[1988]  and  Herbst  [1978],  respectively.  The  species  which  may  be 

vibrationally  or  electronically  excited  are  designated  by  "v"  or 

respectively.  For  reaction  (24),  we  assume  a  10%  yield  of  0(^D).  Reaction 

11  3 

(26)  produces  oxygen  in  the  0(  S),  0(  D)  and  0(  P)  states  with  respective 

branching  ratios  of  0.1,  0.9,  and  1.0  [Zipf,  1970  and  Zipf,  1980].  By 

analogy,  we  assume  that  (27)  yields  atomic  oxygen  with  branching  ratios  of 
0.05,  0.45,  and  0.50. 

C.  Neutral-Neutral  Reactions 

Even  though  this  study  concentrates  on  interactions  between  the 
plasma  and  the  injected  neutrals,  neutral  chemistry  cannot  be  neglected. 
The  neutral  ambient  atmosphere  between  300  and  3000  km  altitude  is  composed 
of  potentially  reactive  species  such  as  0,  N2,  0^,  and  H.  Reactions 

between  atomic  oxygen  and  the  released  gases  are 

0  +  H20  *♦  OH  +  OH  -  0.73  eV  k3A  =  6.0  x  10"15  (38) 

0  +  H20(v)  -t  OH  +  OH  k'3A  =  1  x  10'11 

0  +  CH30H  -*■  H20(v)  +  CH20(v)  +  4.28  eV  k3B  =  8.5  x  10-13  (39) 

All  of  the  reactions  for  this  section  are  given  by  Schofield  [1967]  and 
Kondratiev  [1972]  for  neutral  temperatures  of  1000  K.  Reaction  (38)  is 
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endothermic  and,  consequently,  slow  at  thermospheric  temperatures.  The 

rate  constant  for  reaction  of  f^O  with  atomic  oxygen  is  about  10^  times 

slower  than  the  rate  of  1^0  reacting  with  0+  in  (4).  Above  300  km 

+  -4 

altitude,  the  ratio  0  to  0  concentration  is  10  or  greater.  Therefore, 

thermal  interactions  between  0  and  HjO  are  unimportant. 

For  vibrationally  excited  water,  the  rate  constant  for  (38)  becomes 

k'^A  -  1.0  x  10-11.  A  similar  rate  increase  occurs  for  neutral  f^O  moving 

at  orbital  velocities  relative  to  the  ambient  atmosphere.  At  5  km/s,  the 

center  of  mass  energy  is  2.0  eV  and  the  reaction  cross  section  for  (38)  is 

-7  2 

approximately  1  x  10  cm  [Fraser  et  al.,  1988].  In  this  study,  we  will 
neglect  the  high  energy  reaction  chemistry  and  only  consider  low-speed, 
thermal  reactions. 

Reaction  between  atomic  oxygen  and  water  vapor  will  be  greatly 
enhanced  for  excited  states.  For  example,  metastable  oxygen  reacts  rapidly 
with  H20 

0(XD)  +  H20  +  0H(v)  +  0H(v)  +  1.23  eV  k3c  =  2.0  x  10-10  (40) 

where  OH  can  be  produced  in  the  v=0,l,or  2  vibrational  states. 

Reaction  (39)  is  exothermic  under  all  conditions.  It  is  dominant 
where  k.3B[0]  is  greater  than  or  equal  to  the  quantity  klc(0+].  The 
reaction  products  of  (39)  are  in  vil rational  states  which  will  emit 
infrared  radiation. 

The  following  neutral  reactions  help  determine  the  fate  of  the 
products  from  the  ion-molecule  and  dissociative-recombination  reactions: 

0  +  OH  -»  02*  +  H  +  0.73  eV  k3£)  =  5.0  x  10-11  (41) 

0  +  0H( v)  -»  02*  +  H 


7 


(42) 


OH  +  OH  -»  H20(v)  +  0  +  0.73  eV 
0H(v)  +  0H(v)  -►  H20(v)  +  0 

OH  +  CHa  -»  CH3  +  H20  +  0.75  eV 

0  +  CH3  -»  CH20(v)  +  H  +  3.00  eV 

OH( v)  +  {02  or  N2)  -»  OH(v-l)  +  (02  or  N2) 


k3£  =  2.0  x  10"12 
k'3E  •  3.0  x  10-11 

K3F  =  3.0  x  10"11  (43) 

k3G  «  3.0  x  10"11  (44) 

k3H  =  v  x  8  x  IQ"13  (45) 


Both  (41)  and  (42)  can  produce  significant  reductions  in  the  OH 
concentrations.  The  rate  constants  for  vibrational  relaxation  of  OH  in 
reactions  (41),  (42)  and  (45)  are  provided  by  Llewellyn  et  al.  [1978). 

No  other  neutral  reactions  are  used  in  the  model  because  either  the 
reactions  are  slow  or,  for  H  and  02,  the  background  concentrations  are  low. 


III.  EMISSIONS  FROM  EXCITED  SPECIES 

For  the  water  release,  the  most  likely  sources  of  enhanced  emissions 

are  H20+  [reactions  (4)  or  (5)],  OH  [reactions  (8),  (12),  (13),  (24),  or 

2  2 

(25)],  and  H20  [reaction  (20) J.  Many  bands  of  the  A  A^  ■+  X  electron 
emission  spectrum  of  H20+  occurring  in  the  4000-7000  Angstrom  wavelength 
region  have  been  analyzed  [Lev,  1976).  These  include  bands  which  have  been 
seen  in  the  tails  of  comets  but  were  not  seen  during  the  low-velocity  H20 
releases  (called  Lagopedo)  conducted  over  Hawaii  in  1977. 

Observations  of  the  3064  A  emission  from  the  lowest-energy ,  vibronic 
state  of  OH  were  rep<-  ed  by  Yau  et  al.  [1985)  from  the  H20  release  during 
Waterhole  I.  The  emission  reaction  is  0H(A2E)-K)H(X2Il)+hv(3064  A)  where 
both  the  initial  and  final  states  are  in  the  v  =  0  vibrational  state.  The 
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emission  intensity  was  about  10  Rayleigh.  Only  5%  of  the  OH  in  reaction 

2 

(24)  was  formed  in  the  Z  state  (Anderson  and  Bernhardt,  1978].  The 

2 

important  emission  bands  from  0H(A  Z)  and  their  spontaneous  emission  rates 
are  shown  in  Table  I.  The  initial  and  final  vibrational  levels  are 
indicated  by  (v',  v")  (Crosley  and  Lengel,  1975). 

Table  I 


Vibronic  Transitions  between  0H(A^E+)  v'  and  OH(X^n^)  v" 


(Wavelengths 

in  Angstroms 

and  A  ,  ,, 

v'  ,v" 

"In 
ms  ) 

V" 

\r  f 

0 

1 

2 

3 

4 

V 

0 

3064  , 

3428  , 

1.2x10° 

4.4x10 

1 

2811 

3122  , 

3484  , 

2. 1x10° 

1.3x10 

1.4x10* 

2 

2609  , 

2875  , 

3185 

3550 

4.2X103 

3.3x10 

1.8x10° 

2.2x10 

3 

2444 

2677  , 

2945  . 

3254 

3353 

1.4xl04 

1.7xl03 

S.SxlO^ 

2.6x10^ 

? 

4 

2517 

2753 

3022 

3331 

? 

? 

? 

? 

2 

The  ground  state  of  0H(X  II)  can  be  vibrationally  excited  by  reactions 
(7 ),  (8),  (24),  and  (25).  The  vibrational  and  rotational  bands  for  the 
ground  state  of  OH  extend  from  44745  to  3810  Angstroms  as  shown  in  Table 
II.  The  vibrational  transition  probabilities  have  been  theoretically 
derived  by  Mies  [1974].  We  may  be  able  to  estimate  the  intensities  of  the 
OH  lines  by  comparison  with  the  well  known  reaction  of  ozone  with  atomic 
hydrogen  O^+H-K^+O8*^^  eV  [McDade  and  Llewellyn,  1988]. 


Table  II 

2 

Vibrational  Transitions  in  the  OH(X  IL)  state  (Angstroms) 


..  t 

v"  0 

1 

2 

3 

4 

5 

6 

7 

8 

V 

1 

28016 

2 

14342 

29380 

3 

9791 

15052 

30862 

4 

7461 

10286 

15830 

32495 

5 

6137 

7849 

10684 

16690 

34308 

6 

5254 

6464 

8278 

11285 

17647 

36334 

7 

4627 

5542 

6827 

8758 

11961 

18734 

38668 

8 

4164 

4890 

5865 

7239 

9305 

12743 

20003 

41440 

9 

3810 

4409 

5187 

6234 

7712 

9942 

13662 

21514 

44745 

Transition  Probabilities  for 

0H(X2H 

j)  (s_1) 

v"  0 

1 

2 

3 

4 

5 

6 

7 

8 

V 

1 

20.00 

2 

13.98 

25.23 

3 

0.91 

39.63 

21.35 

4 

0.08 

4.27 

72.02 

13.47 

5 

0.05 

0.39 

10.56 

107.4 

06.79 

6 

0.00 

0.05 

1.27 

20.88 

139.90 

6.03 

7 

0.00 

0.00 

0.18 

2.88 

37.18 

161.0 

14.32 

8 

0.00 

0.00 

0.03 

0.56 

5.64 

60.62 

163.7 

32.41 

9 

0.00 

0.00 

0.00 

0.13 

1.16 

10.81 

89.86 

143.1 

58.12 

The  excitation 

of  atomic 

oxygen 

by  reaction 

(27)  has 

been  not 

been 

verified  by  laboratory  or  space  experiments.  It  is  likely,  however,  that 
the  metastable  0(2D)  and  0(2S)  will  be  produced.  These  states  radiate 


visible  emissions  by  the  following  reactions: 


0(XD) 

-»  0(3P) 

+  hv(6300  A) 

Aj  =  5.1  x  10  3  s  2 

(46) 

0(1D) 

■*  0(3P) 

+  h\>(6363  A) 

A2  =  1.6  x  10”  3  s-1 

(47) 

0(2S) 

-»  0(2D) 

+  h\>(5577  A) 

A3  =  1.34  s'1 

(48) 

Quenching  of  0(2D)  by  collisions  with  molecules  (N^,  C>2,  H^O,  etc.)  is 

represented  by  the  reaction: 

10 


and  by  (40). 

The  H2O  excitation  by  reaction  (25)  has  not  been  experimentally 
verified.  From  energy  considerations,  ve  know  that  this  reaction  may  only 
yield  ground  state  H20(X*A^)  because  the  next  highest  level  is  a  triplet 
state  which  requires  an  energy  greater  than  6.5  eV  for  excitation  [Vang  et 
al.,  1977].  The  H20(X*A^)  vibrational  levels  are  given  by  the  formula 


G(vrv2,v3)  .  oijVj  ♦  UjVj  ♦  »3v3  ♦  XjjVj2  .  X22v22  .  X33*32 


.  X12VjV2  ♦  X13VjV3  ♦  X23v2v3 


where  v^,  and  v^  are  the  vibrational  quantum  numbers  and  the 

coefficient  values  (cm~^)  are 


(Oj  -  3693.8 

«2  -  1614.5 

w3  =  3801.7 

Xn  *  -43.8 

X22  =  -19.5 

X33  =  ”46*3 

X12  =  -20.0 

X13  =  -155.0 

X23  =  -19.8 

The  H2O  spectrum  is 

found  from 

the  difference  in  energy 

levels  by  the 

formula 

X(  Angstroms)  = 

108/[G(v1',v2 

',v3')  -  G(v1",v2",v3")) 

(51) 

The  HjO  emission  bands  and  their  associated  excitation  energy  are  given  in 
Table  III. 
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Table  III 


H„0  Emission  Bands 


(vl*v2,v3) 


Wavelength  (A)  v2  v3 

62696  010 

31732  020 

27385  100 

26625  001 

18755  011 

14548  021 

13790  1  0  1 

11355  111 

9422  201 

9064  003 

8230  2  1  1 

7959  013 

7230  3  0  1 

6984  1  0  3 

6516  3  1  1 

6316  1  1  3 

5945  321 

5918  4  0  1 

5716  203 


■*  (0,0,0)  Transitions 


Energy  (eV)  Spectral  Range  (A) 

0.198 

0.391 

0.453 

0.466 

0.661 

0.852 

0.899 

1.092 

1.316 

1.367 


1.507 

7993 

to 

8448 

1.558 

7846 

to 

8093 

1  ^15 

7059 

to 

7408 

1./75 

6845 

to 

7128 

1.903 

6408 

to 

6626 

1.963 

6275 

to 

6375 

2.086 

5828 

to 

6019 

2.095 

5830 

to 

5999 

2.169 

5665 

to 

5769 

The  transition  probabilities  for  the  ground  vibrational  states  of  1^0 
are  given  in  Table  IV. 

The  emissions  form  the  methyl  alcohol  vapors  will  most  likely  come 
from  chemical  excitation  of  formaldehyde  (CH20)  by  reactions  (18),  (30), 
and  (31).  The  transition  probabilities  for  the  six  fundamental  vibrational 
states  of  CH20  are  given  in  Table  IV. 


12 


Table  IV 

Transition  Probabilities  of  Vibrationally  Excited  Molecules 


Molecule 

Vibrational 

Fundamental 

Transition 

State 

Wavelength 

(A) 

Probability 

A  ,  „  (s-1 

v' t  v" 

h2o 

3(B1) 

26625 

76 

2(Aj) 

62696 

16 

ch2o 

6(B2) 

79554 

2.0 

5(B2) 

35125 

130 

4(Bj) 

84746 

1.7 

3(Aj) 

66357 

2.6 

2(A:) 

57208 

22 

1(AX) 

35971 

124 

IV.  CONTINUITY  EQUATIONS 

The  ionospheric  perturbations  and  airglow  enhancements  can  be 
estimated  by  solving  coupled  continuity  equations  for  the  ions  and 
neutrals.  The  diffusion  equation  for  the  neutrals  is 

=  V-(DVn)  +  P  -  L  (52) 

where  n  is  the  neutral  gas  concentration,  D  is  the  diffusion  coefficient,  P 
and  L  are  the  production  and  loss  terms  representing  the  chemical 
reactions.  The  depletion  of  the  injected  gas  by  diffusion  is  much  more 
rapid  than  the  depletion  by  chemical  reactions  [Bernhardt,  1979].  For  that 
case,  P  and  L  can  be  neglected  and  the  analytic  solution  to  (52)  can  be 
written  as 

N„  o 

n  =  - rjj  exp(-r  /4Dt)  (53) 

(4nDt) 
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where  N  is  the  number  of  Hn0  or  CH-OH  molecules  released  and  r  is  the 

o  15 

distance  from  the  center  of  a  point  release. 

Initially,  our  model  will  neglect  diffusion  of  the  reaction  products 
from  the  chemical  reactions.  This  typically  yields  a  factor  of  two 
overestimate  in  the  airglow  intensities  [Bernhardt,  1984].  The  coupled  set 
of  rate  equations  for  the  densities  are  of  the  form 

3n . 

=  **2 ,  n j )  for  i=l,2,  ...,J  (54) 

where  n^  is  the  density  of  species  number  "i"  and  is  the  production  and 
loss  function  representing  chemical  interactions  with  all  of  the  densities, 
and  J  is  the  total  number  of  species  to  be  calculated. 

The  intensity  is  calculated  with  a  line  of  sight  integral  of  the 
volume  emission  rate.  In  Rayleighs,  the  intensity  is  found  from 

CO 

I  *  10“6  f  A  n  dr  (55) 

io 

where  A  is  the  Einstein  coefficient  of  the  desired  transition  of  excited 
species  n. 

The  time  dependent  equations  are  solved  using  the  program  CHEMGLOU 
developed  at  NRL.  Equation  (54)  is  solved  in  steady  state  for  reactions 

_3 

with  fast  time  constants  (~  10  sec  or  less).  The  inputs  to  the  program 
are  background  neutral  and  ion  densities,  background  temperatures,  and 
number  of  and  CH^OH  molecules  released.  The  outputs  are  temporal 

variations  in  densities  of  ions  and  excited  species,  and  in  airglow 
emission  intensities.  The  program  is  written  in  Fortran  for  use  on  a  VAX 
computer. 
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V. 


SIMULATION  RESULTS 


The  theory  developed  in  the  previous  sections  is  use  to  model  the 

disturbances  produced  by  the  release  of  HjO,  and  CH^OH  into  the  ionosphere. 

♦*  +  6  3  3 

The  ambient  0+  and  H+  densities  are  10°  and  lO-3  cm-  ,  respectively.  The 

26 

simulation  considers  the  release  of  10  molecules  (a  fev  kilograms)  of 

each  species  at  300  km  altitude.  The  neutral  composition  at  that  altitude 

is  given  in  Table  V.  The  computed  diffusion  coefficients  are  3.4  x  10*^, 
10  2 

and  1.0  x  10  cm  /s  for  water  vapor  and  methyl  alcohol,  respectively.  The 
release  is  assumed  to  occur  at  night  so  that  solar  illumination  does  not 
affect  the  chemistry. 


Table  V 

Temperature  and  Concentrations  of 
Neutral  Species  at  300  km  Altitude 


T 

[0] 

to2) 

[No] 

[He] 

[H] 


990 

K 

6.2 

X 

co 

o 

* — 4 

-3 

cm 

4.7 

X 

106 

-3 

cm 

1.3 

X 

>-* 

o 

00 

-3 

cm 

3.7 

X 

106 

-3 

cm 

1.7 

X 

105 

-3 

cm 

The  major  uncertainty  is  the  branching  ratios  for  production  of 

excited  vibrational  states.  For  this  simulation,  we  have  assumed  that  all 

allowed  vibrational  levels  are  equally  populated.  For  example,  all  10 

2 

vibrational  states  of  0H(X  IL  ,  v=0, 1 , 2, . . . , 9)  in  reaction  (19)  are  taken  to 
be  generated  with  probability  equal  to  1/10.  These  states  will  also  be 
generated  by  cascading  from  higher  energy  electronic  or  vibrational  states. 
The  branching  ratios  are  parameters  which  can  be  adjusted  according  to 
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future  theoretical  and  experimental  considerations.  It  should  also  be 
noted  that  rotational  splitting  of  the  lines  has  been  neglected.  This  can 
readily  be  added  for  a  given  rotational  temperature. 

Figures  1  through  6  show  the  computed  composition  out  to  radius  of  19 
km  from  the  center  of  the  release  after  128  seconds.  The  radial  profiles 
of  the  background  atomic  ions  show  a  depression  within  4  kilometers  of  the 
release  (Figure  1).  The  molecular  ions  (0H+,  ant*  ®30+)  ^rom 

water  interactions  fill  in  the  void  produced  by  depletion  of  the  atomic  (0+ 
and  H+)  ions.  Dissociative  recombination  of  the  molecular  ions  and 

H^0+)  is  the  primary  source  for  the  vibrationally  excited  OH  shown  in 
Figure  2.  The  metastable  states  of  atomic  oxygen  are  present  only  because 
they  have  long  lifetimes.  The  intensities  from  metastables  are  much  lower 
than  the  intensities  from  the  molecular  species.  The  primary  source  for 
vibrationally  excited  water  is  reaction  (39)  of  atomic  oxygen  with  methyl 
alcohol.  Only  10%  of  the  H20(v)  comes  from  the  dissociative  recombination 
reaction  (25)  which  can  be  traced  back  to  the  water  release.  The  final 
product  of  the  H2O  release  are  debris  neutrals  which  are  dominated  by  H2O, 
H,  and  0  (Figure  3).  The  final  abundance  of  OH  is  relatively  low  because 
of  its  high  reactivity. 

Seven  organic  molecular  ions  are  produced  as  a  result  of  interactions 
with  the  methyl  alcohol  vapors  (Figure  4).  Each  one  of  these  is  a 
potential  source  of  a  radiation  through  dissociative  recombination 
reactions.  As  shown  in  Figure  5,  the  two  dominant  infrared  emitters  will 
be  vibrationally  excited  formaldehyde  (CHjO)  and  dimethyl  ether  (CH^OCH^). 
The  primary  sources  for  (^(Kv)  are  ion-molecule  reaction  (18), 
dissociative  recombination  reaction  (30)  and  neutral  reaction  (39).  It  is 
likely  that  any  of  the  eight  organic  species  shown  in  Figure  6  will  also 
contribute  to  the  emission  spectra. 
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A  time  history  of  the  airglov  from  the  vibrational,  electronic  and 
vibronic  states  of  several  species  is  summarized  in  Figure  7.  The  sum  of 
the  intensities  from  all  excited  states  of  each  molecule  plotted  as  a 
function  of  time  after  release.  The  superscript  *  and  the  (v)  designate 
vibronic  and  vibrational  excitation,  respectively.  Most  of  the  emissions 
peak  within  a  few  seconds  after  release.  The  metastable  states  of  oxygen 
produce  emissions  peak  much  later.  0H(v),  H^OCv),  and  CHjCKv)  are  the 
strongest  infrared  emitters.  The  most  strongly  excited  of  the  vibronic 
species  will  be  OH*  which  is  produced  by  reactions  (24)  and  (25).  The 
emissions  from  vibronic  states  of  water  are  relatively  weak  because  of  the 
low  percent  of  ^0*  produced  by  dissociative  recombination  of  H^0+. 

Emissions  from  the  metastable  states  of  oxygen  are  the  strongest 
lines  produced  by  CX>2  releases  in  the  ionosphere  [Bernhardt,  1987].  The 
emissions  from  0(^D)  and  0(*S)  are  relatively  weak  for  H2O  and  CH^OH 
releases  because  there  is  no  direct  process  to  excite  them. 

The  spectra  from  excited  OH  covers  240  nm  to  4400  nm  wavelengths. 
The  near  ultraviolet  spectrum  from  the  0H(A  t+)  is  computed  for  four 
seconds  after  release  to  give  the  intensities  in  Figure  8a.  The 
intensities  of  around  1  kR  are  consistent  with  measurements  made  by  a 
rocket  borne  photometer  at  306.4  nm  during  the  Waterhole  I  experiment  [Yau 
et  al.,  1985]  when  H2O  was  explosively  deposited  into  the  auroral 
ionosphere.  Our  computations  show  that  the  emissions  from  vibrational  OH 
will  be  substantially  stronger. 

The  vibrational  spectrum  for  chemically  induced  Meinel  bands  of  OH 

4 

show  intensities  between  1-R  and  3  x  10  -R.  The  visible  emissions  four 
seconds  after  release  are  relatively  weak  (Figure  8b)  because  the 
vibrational  transitions  with  Av  of  4  or  5  are  not  fed  by  cascading  from 
vibronically  excited  states  and  the  transition  probabilities  are  low.  The 
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strongest  emissions  (Figures  8c  and  d)  result  from  low  lying  vibrational 
levels  from  small  values  of  Av. 

After  128  seconds,  the  intensities  of  the  emissions  lines  has  been 
reduced  but  the  levels  are  still  significant.  The  near  ultraviolet 
intensities  are  40-R  or  lower  (Figure  9a)  and  the  Meinel  band  intensities 
(Figures  9b, c  and  d)  are  reduced  by  a  factor  of  thirty  below  the 
corresponding  values  at  t  =  4  seconds. 

VI.  Conclusions 

The  chemical  interactions  of  a  single  molecular  species  with  the 
plasma  constituents  of  the  upper  atmosphere  can  generate  a  large  number  of 
new,  excited  species.  An  1^0  release  by  itself  yields  vibrational  states 
of  H^0+,  ,  0H+,  OH,  and  H2O  along  with  electronic  states  of  OH  and  0. 
Several  sources  of  energy  are  available  to  excite  the  molecules  including 
chemical  energy  from  reactions,  kinetic  energy  from  the  gas  injection, 
solar  energy  and  energy  from  collisions  by  electrons.  This  study  has 
concentrated  on  ion-molecule  and  neutral-molecule  chemistry  as  the  source 
of  excitation  energy.  With  these  sources,  we  find  that  a  wide  range  of 
light  emissions  can  be  stimulated  by  a  single  H2O  release.  Intensities  of 
1  kR  or  more  can  be  generated  between  240  nm  and  4400  nm  wavelengths. 

The  addition  of  CH^OH  to  the  release  produces  many  new  excited 
species  of  the  form  caHjj°c  where  a,  b,  and  c  are  small  integers.  Enhanced 
vibrational  emissions  from  0H(v)  are  produced  by  the  release  of  either 
species.  Emissions  from  C^O  are  uniquely  identified  with  the  CH^OH 
release  but  the  addition  of  water  enhances  the  formation  of  vibrational 
formaldehyde  via  reaction  (18).  The  example  considered  here  demonstrates 
that  the  interaction  of  gases  in  a  mixture  may  enhance  the  emissions  from 
selected  species  (ie.,  CHjO). 
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For  releases  at  300  km  altitude,  the  airglov  emissions  from  chemical 
reactions  may  have  equal  contributions  from  neutral-neutral  reactions,  ion- 
neutral  reactions,  and  electron-ion  reactions.  The  emission  chemistry  for 
releases  much  below  this  altitude  will  probably  be  dominated  by  neutral- 
neutral  reactions  because  of  the  increase  in  the  background  neutral 
densities.  Releases  above  this  altitude  will  be  affected  primarily  by 
reactions  with  the  ambient  plasma.  At  300  km  or  above,  the  loss  of  excited 
states  is  primarily  by  radiation.  At  much  lower  altitudes,  quenching  of 
excited  states  by  reactions  with  neutrals  may  be  important. 

Finally,  our  test  case  has  led  us  to  postulate  that  the  emissions 
from  molecules  released  into  the  nighttime  F-layer  will  come  from  fragments 
of  the  parent  molecule  or  from  excited  atomic  oxygen.  The  1^0  and  CH^OH 
releases  studied  here  yield  enhance  OH  and  C^O  emissions,  respectively. 
For  similar  species,  the  emissions  will  come  from  radicals  produced  by 
detaching  hydrogen  from  the  released  molecule. 

The  next  logical  step  for  this  study  is  to  make  optical  observations 
of  H2O  and  other  releases  which  occur  in  the  night  ionosphere.  Ground 
based  spectrographs  can  be  employed  to  record  the  wide  range  of  emissions 
which  we  have  predicted.  The  observations  should  be  made  under  different 
conditions  of  solar  illumination,  at  different  altitudes,  and  at  different 
injection  velocities  to  separate  the  effects  of  the  different  energy 
sources  described  above.  Such  experiments  have  been  proposed  by  us  for 
future  Space  Shuttle  missions. 
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dissociative  recombination  of  electrons  with  the  molecular  ions 
found  in  Figure  1.  Vibronic  states  of  OH  and  H«0  are  not  shown. 
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Figure  3.  Enhanced  concentrations  of  neutrals  resulting  from  the  water 
component  of  the  release. 
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Figure  6.  Artificially  created  concentrations  of  neutral  molecules  of  the 
form  CHn  (n=2,3,4),  CfM)  (n=0,l»2,3,4),  and  CH3OCH3. 
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Figure  7.  Variations  in  the  light  intensities  from  vibrational  and 
vibronic  states.  Ion-molecule  and  electron-ion  reactions 
between  the  ionosphere  and  the  released  gases  provide  the 
excitation  energy. 
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Figure  8a.  Computed  intensities  of  near  ultraviolet  emissions  from 

rs  n 

OH(A  i+;v')  ->  OH(X^IL,v")  transitions,  four  seconds  after  the 
release. 
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Figure  8b.  Computed  intensities  of  visible  emissions  from  vibrational 
0H(X  n.;v'-*v")  transitions,  four  seconds  after  the  release. 
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Figure  8c.  Computed  intensities  of  near  infrared  emissions  from  vibrational 
transitions  of  OH,  four  seconds  after  the  release. 
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Figure  8d.  Computed  intensities  of  infrared  emissions  from  vibrational 


transitions  of  OH,  four  seconds  after  the  release. 
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Figure  9a.  Computed  intensities  of  near  ultraviolet  emissions  from 
2  +  2 

0H(A  £  ;v' )  ■*  OH(X  n^v")  transitions,  128  seconds  after  the 
release. 
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Figure  9b.  Computed  intensities  of  visible  emissions  from  vibrational 
2 

OH(X  n,;v'-»v")  transitions,  128  seconds  after  the  release. 
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Figure  9d. 


Computed  intensities  of  infrared  emissions  from  vibrational 
transitions  of  OH,  128  seconds  after  the  release. 
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